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ABSTRACT 


A  continuum  model  of  crack  extension  is  used  to  derive  a  crack  propaga¬ 
tion  law  for  the  case  of  constant  plastic  strain  amplitude  fatigue.  The  rate 
of  crack  growth  is  found  to  be  proportional  to  the  square  root  of  the  crack 
length.  Integration  over  the  total  number  of  cycles  to  failure,  Nf  ,  yields 
an  expression  of  the  form  Nf(A«p)n+1  =  f(Aep)  ,  where  fiZ  is  the  applied 
plastic  strain  range,  and  f  is  a  function  which  varies  rapidly  in  the  re¬ 
gion  of  large  strains,  but  approaches  a  constant  as  Aep  becomes  small.  The 
strain  hardening  coefficient,  n  ,  and  the  fracture  strain  enter  as  material 
constants.  Comparison  with  experimental  data  gives  good  agreement  for 
£  0.01  ,  which  is  consistent  with  the  assumptions  used  in  the  theory. 

A  discussion  is  given  which  interprets  the  well-known  power  law 
Nf  d Sp  «  const,  in  terms  of  crack  propagation. 
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I .  INTRODUCTION 

\ 

The  process  of  fatigue  failure  is  normally  divided  into  crack  nucleation 
and  crack  propagation  stages.  Porsythei'  has  shown  that  the  latter  process 
is  itself  composed  of  at  least  two  sequential  stages  which  he  has  termed  Stage 
I,  slip  (shear)  plane  cracking,  and  Stage  II,  tensile  mode  cracking.  Stage  II 
propagation  takes  place  on  a  plane  oriented  90°  to  the  tensile  stress  and  the 
crack  tip  deformation  leaves  characteristic  striations  on  the  fracture  sur¬ 
face.  Each  striation  has  been  shown  to  correspond  to  a  single  stress  (or 
strain)  cycle.  The  study  of  fracture  surfaces  has  revealed  the  existence  of 
these  striations  (or  ripples)  on  a  wide  variety  of  materials  and  demonstrates 
the  generality  of  the  Stage  II  process . 

The  Stage  II  mechanism  has  been  studied  by  Laird  and  Smith§/  in  aluminum 
fatigued  at  high  stress  levels.  By  sectioning  specimens  strained  to  various 
stages  of  the  stress  cycle,  it  was  possible  to  demonstrate  directly  the  crack 
tip  deformation  which  leads  to  the  formation  cf  striations.  The  crack  is  ex¬ 
tended  and  blunted  on  the  tensile  half -cycle,  then  resharpened  on  the  com¬ 
pressive  half-cycle.  The  latter  process  restores  the  stress -concentration 
factor  necessary  for  repetition  of  the  cycle.  It  Is  this  mechanism  which  has 
suggested  the  model  used  in  this  report. 

The  fraction  of  total  life  spent  in  Stage  II  propagation  increases  as  the 
applied  stress  (or  strain)  amplitude  increases . .5/  Laird  and  Smithi/  have  cor¬ 
related  the  number  of  ripples  (observed  optically)  on  the  fracture  surfaces  of 
A1  and  Ni  with  total  cycles  to  failure  and  find  that  for  lives  less  than  100 
cycles,  at  least  70  per  cent  of  life  is  occupied  by  Stage  II.  For  lives  of 
~  103  cycles,  at  least  half  the  life  consisted  of  Stage  II  propagation.  Re¬ 
cent  studies^/  of  a  variety  of  materials  have  confirmed  the  fact  that  crack 
growth  by  ripple  formation  occupies  greater  than  75  per  cent  of  life  in  high 
strain  -  low  cycle  fatigue. 

Nucleation  of  a  propagating  Stage  II  crack  therefore  requires  something 
less  than  25  per  cent  of  total  life  in  low  cycle  fatigue.  At  the  crack  origin 
the  narrow  spacing  ( ~  0.1  p)  of  the  ripples  and  possible  obliteration  due  to 
"banging"  of  the  crack  surfaces  make  It  difficult  to  place  a  definite  lower 
limit  on  the  nucleation  time.  Nevertheless,  in  light  of  the  work  just  cited, 
it  appears  logical  to  approach  the  problem  of  low-cycle  fatigue  in  terms  of 
Stage  II  crack  propagation  and  a  detailed  analysis  of  the  crack-tip  deforma¬ 
tion  during  any  given  cycle. 


Manuscript  released  by  the  author  December  1964  for  publication  as  a  RTD 
Technical  Report. 
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The  treatment  to  be  given  in  this  report , while  based  on  continuum  theory, 
makes  no  pretense  for  elegance. .  Rather,  it  seeks  to  develop  a  simple  model, 
based  on  experimental  evidence  and  plausible  assumptions,  such  that  the  math¬ 
ematics  remains  tractable  and  the  results  capable  of  experimental  verifica¬ 
tion.  The  case  of  fully  reversed,  constant  plastic  strain  amplitude  was 
chosen  both  for  its  simplicity  and  the  fact  that  a  well-documented  empirical 
law  exists  which  relates  the  number  of  cycles  to  failure  to  the  applied  strain 
amplitude 


II.  CRACK  GROWTH  MODEL 


The  experimental  results  of  Laird  and  Smith£/  on  aluminum  and  of 
McEvily  et  al.— '  on  polyethylene  (which  also  produces  Stage  II  striations ) 
have  been  used  to  construct  the  continuum  model  of  crack  extension  illustrated 
in  Figure  1.  The  crack  is  a  two-dimensional  ellipse  and  considered  to  be 
located  somewhere  in  a  thick  specimen  such  that  conditions  of  plane  strain 
exist  at  the  tip.  In  the  fully  compressed  state,  the  crack  has  semi -major 
axis  equal  to  c  and  semi -minor  axis  equal  to  b0  .  As  the  applied  strain 
increases  to  the  maximum  tensile  state,  plastic  deformation  at  the  crack  tip 
results  in  the  creation  of  new  surface  and  the  crack  assumes  the  shape  with 
semi-major  axis  equal  to  c  +  6c  ( 6c  =  RR 1  )  and  semi -minor  axis  equal  to 
bRl  .  In  the  derivation  to  follow,  it  is  assumed  that  6c  «  c  .  During  the 
compression  half -cycle,  the  crack-tip  strain  reached  during  tension  is  con¬ 
sidered  irreversible  such  that  the  arc  P*R*  becomes  P"R"  in  the  fully 
compressed  state.  In  this  manner  an  extension  equal  to  Ax 2  is  achieved 
during  one  complete  cycle.  The  folding  which  leads  to  ripple  formation  has 
been  neglected  in  the  calculation  of  Ax2  . 


The  assumption  of  strain  irreversibility  at  the  crack  tip  is  an  example 
of  the  instability  of  ductile  materials  under  reversed  loading  which  has  been 
discussed  by  MeClintock.  9/  In  fact,  the  model  presented  here  can  be  consid¬ 
ered  a  special  case  of  MeClintock* s  "growth  by  progressive  deformation"  case, 
for  which  he  has  shown  the  instability  to  arise  from  either  the  Bauschinger 
effect  or  low  strain  hardening  rates. 


Referring  to  Figure  1,  the  root  strain,  erq  >  is  defined  as 
(OP'-OPVOP  ,  where  0  is  chosen  so  that  OR  is  the  root  radius,  rQ  ,  of 
the  compressed  crack.  To  relate  Ax2  to  this  strain,  consider  the  triangles 
OP^'  and  OP"R"  .  These  right  triangles  are  approximately  equal  since  the 
arcs  P'R*  and  P"R"  are  taken  to  be  equal  and  OR*  &  OP".  Thus  RR"  ~  PPf  ; 
and  since  PP*  can  also  be  expressed  as  OP* -OP  =  r0er}_  ,  the  following  rela¬ 
tion  holds, 
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CRACK  GROWTH  MODEL  , 


Figure  l  -  Continuum  Model  for  S^age  II  Fatigue  Crack  Growth 


=  Fo*rl 


(1) 


where  p  is  a  constant  of  the  order  1  which  allows  Equation  1  to  be  written 
as  an  exact  equality. 


Ill .  DERIVATION  OF  CRACK  GROWTH  IAW 


The  framework  of  the  derivation  is  as  follows;  The  root  strain  is  first 
expressed  in  terms  of  the  applied  strain  through  a  strain  concentration  fac¬ 
tor.  The  strain  hardening  coefficient  of  the  material,  n  ,  is  introduced  at 
this  point.  Next,  the  crack  growth  rate,  dc/dN  «  tag  >  is  integrated  between 
the  limits  N0  ,  the  number  of  cycles  to  initiate  a  Stage  II  crack,  and  Nj>  , 
the  number  of  cycles  to  fracture.  Tfie  limiting  values  of  the  crack  length 
are  cQ  ,  the  initial  Stage  II  crack  length,  and  Cf  ,  the  final  crack  length, 
for  which  the  next  quarter-cycle  of  strain  will  bring  ductile  fracture .  The 
fracture  strain  of  the  material,  ,  is  introduced  at  this  point.  Thus,  two 
material  constants,  n  and  “E^,  ,  appear  in  the  final  result. 

The  root  strain  can  be  expressed  in  terms  of  the  changes  in  the 

semi-minor  axis,  =  (bD -bQ )/bQ  ,  provided  6c  «  c  .  The  calculation  re¬ 
lating  this  strain  to  the  macroscopic  applied  strain  is  more  accurately  car¬ 
ried  out  using  logarithmic  strains  since  the  magnitudes  are  finite.  Thus,  we 
define  the  incremental,  logarithmic  strain  as  d"?^  =  db/b  *  l/2(dr/r)  ,  where 
r  is  the  root  radius .  The  total  logarithmic  root  strain  is  therefore 
Crl  =  l/2  ln(rm/r0)  at  the  conclusion  of  the  tensile  half-cycle. 

To  evaluate  the  strain  concentration  factor,  we  recognize,  as  DruekeriS?/ 
has  suggested,  that  this  factor  will  vary  as  the  crack-tip  geometry  changes 
during  the  cycle.  An  incremental  strain  concentration  factor,  k£  ,  is  there¬ 
fore  defined  by  dc^  -  kcdc  ,  where  e  is  the  applied  logarithmic  strain. 
Neuberii/  has  given  an  expression,  valid  in  shear,  which  relates  the  Hookian 
stress  concentration  factor,  k^  ,  to  the  stress  and  strain  concentration  fac¬ 
tors,  and  kc  ,  for  any  arbitrary  nonlinear  stress-strain  law, 


(Ve)l/2  " 


(2) 


Neuber  has  further  stated  that  the  relation  may  be  extended  to  arbitrary 
states  of  stress  without  modification.  Rhee  and  McClintocki^/  have  questioned 
the  validity  of  this  procedure  in  the  tensile  case  for  nonhardening,  ideally 
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plastic  materials.  In  the  present  calculation,  for  want  of  a  better  relation, 
we  shall  use  the  Neuber  expression  as  valid  for  materials  in  tension  with  fi¬ 
nite  strain  hardening.  Assuming  that  the  power  law,  o  =  oQen  ,  holds  for  the 
materials  in  question,  we  find 


1 

n+1 


where  has  been  specified  for  an  elliptical  crack  of  length  2c  and  root 
radius  r  .  The  strain  hardening  coefficient,  n  ,  should  strictly  be  inter¬ 
preted  as  characteristic  of  the  material  in  the  vicinity  of  the  crack  tip  as 
the  tensile  half-cycle  is  applied  to  it.  In  the  sequel,  because  such  values 
are  unknown,  we  shall  take  n  to  be  the  usual  value  determined  from  a  tensile 
stress-strain  curve. 


Equation  3  may  now  be  used  in  the  incremental  expression  for  dc^  to 
find  the  relation  between  the  root  strain  ani  the  applied  strain.  We  have 


1 

n+1 


Integrating  between  rQ  and  rm  on  the  left,  and  0  and  ,  the  total 

plastic  strain  range,  on  the  right,  we  obtain 


1 


n+1 


1 

n+1 


or, 


€rl 


n+1 

2 


In  a 


(4) 


Here  it  has  been  assumed  that  the  applied  strains  are  large  enough  to  neglect 
the  elastic  component,  and  that  the  total  plastic  strain  range  corresponds  to 
opening  the  crack  from  its  fully  compressed  state  to  its  maximum  opening. 
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We  next  assume  that  at  the  conclusion  of  each  strain  cycle,  the  crack 
is  returned  in  its  fully  compressed  state  to  a  root  radius,  rQ  .  Therefore, 
Equation  4  represents  the  maximum  root  strain  for  any  crack  length  c  .  Con¬ 
verting  Equation  4  to  conventional  strains  gives, 


n+1 


and  thus  the  crack  growth  rate  is 


dc/dN  =  Br0  ya  2  -  1 J  .  (5) 

Equation  5  predicts  that  the  rate  of  \rack  growth  in  Stage  II  is  approximately 
proportional  to  the  square  root  of  the  crack  length  under  the  conditions  of 
constant  plastic  strain  amplitude.  This  dependence  on  crack  length  is  quite 
sensitive  both  to  the  way  in  which  the  incremental  root  strain  is  de¬ 

fined  and  to  the  assumption  about  rQ  .  If,  for  example,  the  root  strain  is 
taken  simply  as  dr/r  ,  then  dc/dN  will  depend  linearly  on  c  .  Moreover, 
if  rQ  actually  increases  as  the  creek  lengthens,  then  Equation  5  is  altered 
so  that  dc/dN  increases  more  rapidly  vith  c  .  For  the  balance  of  this  paper 
the  simple  model  which  led  to  Equation  5  will  be  retained. 


Equation  5  can  be  integrated  between  suitable  limits  to  find  the  number 
of  cycles  spent  in  Stage  II  crack  growth.  The  integration  can  be  expressed 
in  terms  of  a  to  give, 


(o-l  )ndcy 


n+1 

(“2  -1) 


=  8 


2n+3 

(n+"un+2 


(Nrl/4) 


dN 


(6) 


Here  N0  is  the  number  of  cycles  to  initiate  a  Stage  II  crack  of  length  c0  , 
and  Nc  the  total  number  of  cycles  to  failure.  an  and  are  obtained  by 

substituting  the  values  for  cQ  and  cp  into  the  expression  for  a  .  The 
final  crack  length  is  achieved  one  quarter-cycle  before  complete  failure 

by  ductile  fracture.  The  integral  on  the  left  in  Equation  6,  designated  Dy 
I  in  the  sequel,  cannot  be  integrated  in  closed  form  except  for  the  limiting 
cases,  n  =  0,  1  . 


To  choose  a  reasonable  value  for  Cp  ,  we  define  a  critical  root  frac¬ 
ture  strain  such  that  when  c  =  cf  ,  ?rl  =  ?rf  =  Hli  In  af  .  A  limiting 
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condition  is  given  when  the  applied  strain  range  &Cp  is  equal  to  the  frac¬ 
ture  strain  ef  ,  and  cf  is  essentially  the  void  size.  For  convenience ,  we 
take  the  void  size  to  be  of  the  order  of  cG  .  One  then  obtains 


cr  = 


The  limits  aQ  and  Op  are  given  by 


Q'f  « 


o  “ 


(?) 


(8) 


The  only  parameter  to  enter  the  evaluation  of  the  integral ,  I  ,  is  the  ratio 
c0/rQ  ;  the  values  of  n  and  ep  depend  on  the  material  under  test,  and 
is  given  by  the  experimental  conditions. 

Because  we  have  taken  the  macroscopic  fracture  strain  ip  as  the  cri¬ 
terion  for  determining  cp  ,  it  should  be  noted  that  the  larger  the  parameter 
c0/r0  )  the  larger  the  required  root  strain  to  achieve  fracture,  *erp  .  If 
€rf  vere  known,  the  best  value  for  cQ/r0  could  then  be  chosen  without  dif¬ 
ficulty.  As  it  is,  a  choice  for  cG/ro  must  he  tf^de  from  comparison  of 
Equation  6  with  experimental  life  data  in  the  low-cycle  range.  Even  so,  the 
derivation  of  Equation  6  with  only  one  adjustable  parameter  (8  enters  only 
as  a  constant  multiplier)  can  be  considered  fortunate. 

Equation  6  can  be  rewritten  in  the  following  form, 


n+2 

(Nf  -  1/4  -  N0)(A?  )n+1  =  iHiil -  I  (9) 

2n+3p 

If  we  take  N0  to  be  small  compared  with  Np  in  the  high  strain  region,  an 
assumption  consistent  with  the  trend  of  recent  experimental  results,  we  have 


Nf(Aep)n+i  =  0j+.l)n+2  I 
2n+3B 


(10) 
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an  equation  which  bears  some  resemblance  to  the  Manson-Coffin  power  law 
relational!/  for  low-cycle  fatigue.  Weiss!5/  has  derived  a  similar  relation 
starting  from  somewhat  different  assumptions,  but  with  a  constant  term  on 
the  right-hand  side.  On  the  other  hand,  the  right-hand  side  of  Equation  10 
is  a  function  of  Aip  through  the  integral  I  ,  which  approaches  a  limiting 
value  as  A€p - >0  . 

In  the  following  section.  Equation  10  is  compared  with  the  experimental 
results  obtained  by  a  number  of  workers  for  constant  plastic  strain,  low-cycle 
fatigue  tests.  In  all  cases  the  factor  0  was  taken  equal  to  1,  and  the 
integral,  I  ,  was  evaluated  by  Simpson's  Rule  using  an  IH4  1620  automatic 
computer. 


IV.  COMPARISON  WITH  EXPERIMENT 


In  Figure  2,  Aip  versus  curves  for  nickel*  and  three  different 

values  of  the  parameter  c0/r0  are  compared  with  total  life  data  for  Ni  taken 
from  the  work  of  Coffin  and  Tavernelli.ll/  In  these  and  subsequent  curves, 
one-quarter  cycle  has  been  added  to  the  value  of  Nf  calculated  from  Equation 
10  so  that  the  theoretical  and  experimental  data  are  effectively  normalized 
at  the  fracture  strain.  It  is  clear  that  the  value  ^q/^o  =  provides  the 
best  agreement  with  experiment  in  the  region  Ag^  >0.1  where  the  assumptions 
leading  to  Equation  10  should  be  most  applicable'  This  ralue  for  c0/rQ  is 
also  more  satisfying  from  a  physical  standpoint  than  either  of  the  other 
values.  Further  consequences  of  this  choice  for  cQ/r0  will  be  discussed 
at  the  end  of  this  section. 

The  sensitivity  of  the  computed  life  values  to  various  choices  for  n 
and  ef  ,  is  shown  in  Figures  3  and  4  for  nickel.  The  fact  that  n  -  0*4  is 
the  accepted  value  is  quite  clear.  The  curve  for  =  2.0  in  Figure  4  would 
be  the  result  if  A?^  at  one-fourth  cycle  were  chosen  as  instead  of 

,  as  has  been  suggested  by  Willner  and  McClintock.i^/  Clearly  this  choice 
would  not  produce  better  agreement  with  the  experimental  data  in  the  present 
case. 


ef 


*  Sources  of  the  material  constants  n  and 


are  given  in  Table  1. 


TABLE  1 


MATERIAL  CONSTANTS 


Metal 

n 

If 

1100  A1 

0.5^/ 

1.71 

Ni 

0.4^ 

1.10 

2024  A1 

0.2^/ 

0.4 

0FHC  Cu 

0.54^./ 

1.4 

a/  A.  M.  Willner  and  F.  A.  McClintock,  Paper  No.  61-WA-273 
Annual  Meeting  of  ASME  (1961). 

b/  G.  W.  Gel  1  and  N.  L.  Carwile,  NBS  Circular  520,  67  (1952). 
c/  J.  R.  Low  and  F.  Garfalo,  Proc .  Soc .  Exptl.  Stress  Anal. 

4,  16  (1947). 

d/  J.  R.  Low,  "Properties ..of  Metals  in  Materials  Engineering, 
ASTM  ( 1949 ) . 

e/  J.  F.  Tavernelli  and  L.  F.  Coffin,  Jr.,  Trans .  of  ASME 
84  D,  533  (1962). 
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In  Figure  5,  the  results  of  Equation  10  are  compared  with  the  experi¬ 
mental  data  for  1100  A1  from  two  sources  .—til/  Again,  the  choice  of  c0/rQ  = 
10  is  consistent  with  the  data.  Figure  6  compares  the  theory  with  the  data 
for  OPHC  copper,  and  again  shows  the  sensitivity  of  the  theory  to  the  entice 
of  n  (accepted  value,  n  =  0.5).  Figure  7  compares  the  theory  with  data  for 
2024  Al.  This  was  the  only  case  in  which  a  choice  of  cQ/r0  higher  than  10 
was  indicated.  The  low  ductility  of  this  alloy  compared  to  the  other  three 
materials  may  account  for  this  fact. 

An  estimate  of  the  order  of  magnitude  of  c0  and  rQ  can  be  obtained 
from  experimental  observations  on  crack  growth  rates  in  Stage  II.  If  Equation 
5  is  rewritten  in  nonditaensional  form,  (dc/dN)/rQ  may  be  plotted  against  the 
ratio  c/c0  for  various  choices  of  c0/r0  .  Such  a  plot  for  n  =  0.4  and 
Aep  =  0.02  is  shown  in  Figure  8.  If  is  chosen  as  1.5,  the  ratio  Cj/c0 

is  fixed,  and  the  upper  limit  to  (dc/dN)/rQ  can  be  found  for  the  three 
cases  shown.  McEvily  et  al.Jj/  have  found  for  OFHC  Cu(n  =  0.5,  =  1.4) 

that  the  ripple  spacing  Just  prior  to  fracture  is  of  the  order  of  10“  ^  cm. 
when  Ae^  =0.02  .  This  information,  coupled  with  the  upper  limit  for 
(dc/dN)/rc  can  be  used  to  determine  cc  and  r0  as  shown  in  Figure  8.  The 
value  cQ  =  10“ 3  cm.  when  cQ/r0  =  10  is  reasonable  in  the  light  of  obser¬ 
vations  by  Forsythe3/  and  Grosskreutzi3/  at  much  lower  strain  amplitudes  where 
the  transition  from  Stage  I  to  Stage  II  cracking  in  aluminum  occurs  at  crack 
lengths  of  the  order  of  10  p,.  The  value  for  rQ  when  c0/rQ  -  10  seems  a 

bit  large  on  first  glance,  but  if  one  recalls  that  the  original  model  is  an 

idealized  elliptical  crack,  then  rQ  =  10~4  cm.  is  not  unreasonable. 

Finally,  the  ratio  of  final  to  initial  crack  growth  rate  in  Figure  8  is 

38,  which  compares  with  17  for  the  ratio  of  largest  to  smallest  ripple  spac¬ 

ing  found  by  McEvily  et  al.3/  in  Cu.  The  theory  thus  predicts  that  somewhat 
smaller  striations  should  exist  than  those  actually  resolved  by  electron 
microscopy  of  iracture  surface  replicas. 
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V  .  ^ 


V.  DISCUSSION 


The  good  agreement  between  Equation  10  and  the  experimental  data  over 
several  decades  of  life  lends  good  support  to  the  existing  experimental  data 
which  show  that  low-cycle,  high-strain  fatigue  can  be  viewed  mainly  as  Stage 
II  crack  propagation.  The  expression  first  proposed  by  Manson  and  Coffin,^;?/ 
is  therefore  a  crack  propagation  law.  The  cracks  involved  are  never  very 
large,  ranging  from  about  10  jjl  at  the  outset  to  the  order  of  millimeters  near 
the  end  of  life. 

The  theory  which  has  been  described  suggests  that  the  number  of  cycles 
to  initiate  a  Stage  II  crack,  N0  ,  is  indeed  quite  small  compared  to  Np 
until  the  strain  range,  Aep  ,  falls  below  ~0. 01,  or  Np  >  5  x  10^  cycles. 

It  is  below  this  value  of  Alp  that  the  theoretical  curves  in  Figures  2-7 
begin  to  fall  away  from  the  experimental  points.  It  is  surprising  that  the 
agreement  should  hold  out  to  these  relatively  long  lives.  laird  and  Smithf/ 
have  observed  that  only  30  per  cent  of  the  fatigue  life  is  taken  up  with 
ripple-propagation  for  Np  ~  5  x  10 ^  cycles.  However,  their  observations 
were  limited  to  optical  magnifications,  and  it  is  quite  possible,  in  view  of 
later  work, 2./  that  a  much  larger  percentage  would  have  been  observed  at 
electron-optical  levels.  Another  possibility  is  that  the  simple  model  which 
leads  to  Equation  5  is,  in  fact,  not  accurate  enough,  and  that  the  crack 
growth  rate  takes  an  altered  form.  It  has  already  been  suggested  above  that 
some  elementary  changes  in  the  definitions  of  d”erp  and  r0  lead  to  a  linear, 
or  higher  power,  dependence  of  dc/dN  on  c  .  In  this  case,  computed  values 
of  (Np-NQ)  would  fall  away  from  the  actual  life  curves  sooner  than  in  the 
present  case.  There  is  some  recent  evidence  for  a  linear  relation  between 
dc/dN  and  c  ,5/  which  suggests  that  N0  becomes  appreciable  at  larger  values 
of  Aip  than  those  predicted  here. 

The  sensitivity  of  the  theory  to  variations  in  n  and  7p  (Figures  3, 

4  and  6)  argue  that  the  continuum  model  on  which  the  calculation  is  based 
is  realistic.  To  be  sure,  the  model  has  a  "built-in"  inverse  dependence  be¬ 
tween  Np  and  Aep  .  However,  with  only  one  adjustable  parameter,  c0/rQ  , 
it  is  reasonable  to  conclude  that  the  agreement  between  Equation  10  and  the 
experimental  data  is  more  than  fortuitous. 

It  is  also  clear  why  many  materials  behave  so  nearly  alike  under  the 
same  constant  plastic  strain  conditions.  Values  for  n  are  nearly  all  in 
the  range  0.2  -  0.5,  and  most  values  for  cp  are  of  the  order  of  unity. 

Thus,  the  rate  of  ciack  propagation  for  a  given  ATp  and  c  will  not  vary 
much  over  20  per  cent  over  the  range  of  n  . 
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The  theory  here  developed  is  restricted  to  axial  loading  tests  in  which 
a  single  crack  dominates  most  of  the  te^t  life.  It  has  been  shown,  however, 
that  materials  obey  a  power  law  relation  when  fatigued  in  high-strain 
torsion. Microscopic  observation  at  somewhat  lower  strains^?/  has  shown 
that  torsional  fatigue  produces  a  multitude  Stage  I  cracks  at  the  surface 
before  Stage  II  propagation  begins.  A  dilemma  is  posed,  therefore,  in  the 
interpretation  of  high-strain  torsion  tests.  It  would  be  extremely  inter¬ 
esting  to  study  the  fracture  surfaces  of  such  specimens  to  determine  the 
fraction  of  life  spent  in  Stage  II  propagation. 

Finally,  it  should  be  pointed  out  that  the  theory  predicts  (Equation  15) 
that  near  the  end  of  life,  the  ripple  spacing  for  a  given  material  will  be 
independent  of  the  value  of  A?p  .  It  would  be  interesting  to  explore  this 
possibility  experimentally. 
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